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We demonstrate significant modification of the electron-phonon energy loss rate in a many- valley semiconduc- 
tor system due to lattice mismatch induced strain. We show that the thermal conductance from the electron 
system to the phonon bath in strained n+Si, at phonon temperatures between 200 mK and 450 mK, is more 
than an order of magnitude lower than that for a similar unstrained sample. 



The interaction between electrons and phonons is one 
of the most important parameters of a semiconductor sys- 
tem, because it dictates the room temperature electron 
mobility and electron-phonon (e-ph) energy loss. 'I' In the 
low temperature (sub-1 K) regime e-ph coupling typically 
dominates the total thermal coupling from electrons to 
the environment. Therefore, e-ph coupling is one of the 
key parameters of interest in low temperature bolometer 
and microcooler applications.^ In single- valley semicon- 
ductors the electron-phonon interaction can be typically 
described by deformation potential or piezoelectric cou- 
pling constants, which means that the dependence of the 
e-ph matrix elements on the electronic variables, such as 
momentum, can be typically ignored.'^ The situation is 
quite different in metals where momentum dependency 
must be included explicitly. 

Simply setting the e-ph matrix elements to a constant 
is not the whole story even for deformation potential cou- 
pling in semiconductors. Certain lattice perturbations 
can lift the valley degeneracy so that the valley degree of 
freedom plays an important role in low temperature e-ph 
energy loss.SEl Mean- field theory predicts that the e-ph 
energy loss rate is strongly enhanced in many- valley semi- 
conductors in comparison to single- valley ones due to lack 
of screening.^ This suggests that if one is able to depop- 
ulate asymmetric valleys of a many- valley semiconductor 
then a strong reduction in e-ph energy loss rate should 
be observed. In this letter, we will experimentally test 
this prediction by using strain to depopulate the four in- 
plane valleys of epitaxial n+Si films [see Fig.jlja)]. Sub-1 
K e-ph energy loss experiments of strained samples are 
compared to experimental data of Ref. 4 and our non- 
strained control sample. We find that qualitatively the 
experiments agree with theory and e-ph thermal conduc- 
tance is indeed reduced in the strained sample. The effect 
is, however, smaller than ideally expected. 

Strained silicon epi-layers were produced using a RP- 
CVD epitaxial growth system. Initially, a linearly graded 
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Sii-a^Gca; alloy layer of thickness 2 /im is grown, with the 
germanium content (x) graded up to 0.2. This is followed 
by a constant composition layer of 500 nm Sio.8Geo.2. 
Finally, an epitaxial silicon layer of 30 nm is grown, which 
conforms to the lattice spacing of the alloy layer and is 
thus tensile strained. This silicon layer has a phosphorus 
dopant concentration of 4 x 10^^ cm~^. As a control, a 
similar doped silicon layer was grown on silicon without 
a SiGe virtual substrate. The degree of relaxation of the 
SiGe buffer layer and strain in the surface silicon layer 
were determined by X-ray diffraction. We found that 
the strain in the silicon layer is approximately 0.95% and 
hence the band splitting between the 2-fold and 4-fold 
conduction bands should be about 150 meA^, which is 
sufficient that, at the targeted carrier concentration, only 
the 2-fold bands are populated. 

The geometry of the devices fabricated on these wafers 
is presented in Fig. [ijb). Rectangular mesa pillars of 
height 100 nm (active layer 30 nm), were etched using a 
CF4:02 (10:1) plasma to define the geometry of the de- 
generately doped material and isolate the mesas from the 
surrounding area. After native oxide removal in HF, the 
aluminum leads were then sputter deposited to a thick- 
ness of 200 nm and patterned using a phosphoric acid 
based etch solution. 

The degenerately doped silicon region is connected to 
the aluminium (superconductor) leads through Schottky 
barriers.^ These kind of superconductor-semiconductor 
junctions exhibit similar characteristics as the well known 
superconductor-insulator-normal metal junctions.^ The 
ratio of the subgap resistance to the asymptotic resis- 
tance, a common figure of merit for these kinds of junc- 
tions, is 200 for the strained sample and 40 for the con- 
trol. The temperature dependency of the I-V charac- 
teristics of these junctions is used for thermometry in 
our experiment. We bias the two junctions in the mid- 
dle of the bar with a constant current from a floating 
high impedance current source and monitor the induced 
voltage, which depends on the electronic temperature 
of the semiconductor island.^ The bias current through 
the thermometer junctions is selected so that it has a 
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FIG. 1. (Color online) (a) Schematic illustration of the con- 
stant energy ellipsoid of Si conduction band valleys and how 
the strain should affect these. The four in-plane valleys have 
their energy raised and are depopulated due to strain, (b) 
Simplified geometry of the sample and a sketch of the mea- 
surement setup, Al in dark and Si in light color. The center 
n^ Si mesa is 205 /xm long and 5 /im wide. We voltage bias 
the wire from the ends and measure the induced heating cur- 
rent /. Simultaneously the thermometer junctions are biased 
with a constant current and the voltage V is monitored. This 
voltage gives the electron temperature Tg. 



TABLE I. Sample parameters. Junction resistance and Si 
sheet resistance measured at 50 mK, carrier densities and mo- 
bilities from Hall measurements at 10 K. 



Sample 


Strained 


Control 


Si-Al junction resistance (kQ/xm^) 


108 


12 


Si sheet resistance (Q) 


569 


383 


Carrier density (10"*^^ cm"*^) 


2.7 


3.1 


Mobility (cm^V) 


155 


192 


Mean free path (nm) 


5.3 


5.0 



negligible heating effect. We calibrate these thermome- 
ters against the RuO thermometer of our dilution cryo- 
stat. These junction thermometers, however, saturate at 
around 200 mK, because of the relatively low subgap re- 
sistances. Hence, all the data we show are above this 
temperature. Sample parameters are presented in Ta- 
ble ID 

As shown in Fig. [ijb), we induce heating current / 
through the semiconductor island with voltage bias. The 
total heating power P is calculated from P = Rsil^^ 
where Rsi is the resistance of the n+Si mesa (measured 
in a 4-probe setup). We neglect the heating effect from 
the Al-Si junctions at the ends, assuming that because 
of the high resistivity of the n+Si, and the distance of 
the thermometer from the junctions, the total electron- 
phonon coupling is much stronger than the heat conduc- 
tivity through the wire. Hence, the junction heating is ef- 
fectively "screened" by the electron-phonon coupling. We 
checked this assumption with finite-element-method sim- 
ulations taking electronic thermal conductivity from the 
Wiedemann-Franz law and using our measured electron- 
phonon couplings (see below). The effect of junction 
heating on the electronic temperature near the electron 
thermometer was less than 1 % in all cases. 
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FIG. 2. (Color online) (a) Measured electronic tempera- 
ture of the control and strained sample as a function of the 
heating power applied at different bath temperatures, (b) 
The electron-phonon thermal conductance Ge-ph, determined 
from the derivative of data in (a). Circles showing the lowest 
heating power data points, ie, points where Te is very close to 
Tph . Solid line from Ref . (4] Dashed line calculated OFa / dTe 
and dash-dotted and dotted lines calculated (dFs/dTe) x 300 
in 3 and 2 dimensions, respectively. See text for further de- 
tails. 

The main results of the paper are presented in Fig. [2j 
Figure [2| a) shows the electron temperature Tg measured 
with the thermometer junctions as a function of applied 
heating power at different bath temperatures. One can 
readily see that with the same power the strained sample 
is heated to much higher temperatures than the control. 
This indicates that the e-ph coupling has indeed been 
reduced in the strained sample as compared to the un- 
strained control. 

To obtain more quantitative results we extract the e- 
ph thermal conductance Ge-ph = dP/dT^. This is plot- 
ted in Fig. |2|b) as small dots. The difference between 
conductances in the control and in the strained sample 
is about a factor of 20 at the lowest temperatures and 
a factor of 50 at 400 mK. We have assumed that the 
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phonon temperature T^^ equals the bath temperature. 
As this assumption might break down when the tem- 
perature difference between the electrons and phonons 
becomes large, we plot, as larger circles, the lowest heat- 
ing power data points (P ~ 0), where this assumption 
should be most reliable. We also show as a reference the 
result from a comparable unstrained sample from Ref. 4 , 
Note, however, that in Ref. 4 the measurements extended 
only up to 350 mK, the samples had a dopant concen- 
tration of 3.5 X lO^^cm"^ and they were fabricated on 
silicon-on-insulator wafer. 

In order to compare these results with the theory, we 
assume that the electron and phonon sub-systems are 
each in separate internal equilibrium at temperature Tg 
and Tp^, respectively. The total energy loss rate (or 
power loss) P, which is equivalent to the heating power, 
can then be quite generally described by a symmetric 
formula 

P = V(F(Te) - F{T„,)), (1) 

where V is the volume of the electron gas. Hence 
the electron-phonon thermal conductance is given by 
Ge-ph = dF/dTe. According to the mean- field analysis 
of Ref. 5 , in semiconductors, the energy loss rate function 
F can be expressed as a sum of two components 

F{T)= Fs{T)^Fa{T). (2) 

In single-valley systems only the symmetric part Fs ex- 
ists. Certain lattice perturbations lift the valley degen- 
eracy in many-valley systems and this gives rise to the 
asymmetric part Fa- In the diffusive limit, where the 
electron mean free path (/g) is smaller than the inverse 
thermal phonon wave number (qt = ) , we can write^ 



Fs{T) = K,,,, J: V, 4^^^^^ (3a) 

^ A 

Fa{T) = K,kZ ^ ^A (Si) 4- (3b) 

^ A 

Here Km = 21!^% f where p, and vp are the den- 
sity of states (for electrons), mass density of the crys- 
tal and Fermi velocity, respectively. Numerical factor 
Bm = / dxx^/{e^ — 1), the screening wave number (in- 
verse of the screening length) k, = 2 (e^z^/4£5)"^^*^^ 
and Eb is the background dielectric constant. The di- 
mensionality of the electron system is defined through 
the parameter d = 2,3. Note, however, that Fa does 
not depend explicitly on d. 7 and 7 denote the to- 
tal scattering rate and total intervalley scattering rate, 
respectively {1^ = vp/j). Effective deformation poten- 
tial constants and phonon velocity vx depend 
on mode index A. Here X = L denotes longitudinal 
and X = S denotes transverse mode. If we assume 
symmetry of the Si conduction band valleys and bulk 
phonon system with isotropic acoustic modes we have 



=E\^^El) and (S^) = ^S^(f S^) for A = L(^), 
where fS^S^ + 22^ + . Here Srf(^) is the di- 

latational (uniaxial) deformation potential constant.^' We 
have assumed valley degeneracy of two (six) for ^Sl^^^^^y 

Note that in Eqs. ([3| P5 ex whereas Fa 

does not exhibit any dependency on the screening wave 
number This reflects the quite general principle ac- 
cording to which Fs is strongly screened (at low temper- 
atures) and Fa is unscreened.^ Lack of screening in Fa 
follows physically from the many- valley picture, where 
a phonon can create different perturbation to different 
valleys. Therefore, at low temperatures we should al- 
ways have Fa ^ Ps- in a many- valley system and hence 
P Fa- If valleys with different symmetry can be de- 
populated then Fa = ^ and we should observe a large 
reduction in P, because now F = Fs. 

Consequently we would expect the control sample to 
have Ge-ph — dFA/dT^^ whereas in the strained sample 
we should have Ge-ph = dFs/dTe. Figure [2|b) shows 
these theory predictions as lines. All fixed parameters of 
the electron system(s) are obtained from standard single 
electron formulas. For the other parameters we use typ- 
ical values for Si: p = 2.33 x 10^ kgiR-^, Vt{l) = 4700 
(9200) m/s and Srf(^) = -11.7(9.0) eV.^ Fa has one fit- 
ting parameter 7/7, to which we get a value 1.27 from a 
fit to the control sample data. The fit was performed in 
the temperature range 200 — 270 mK. In this range the 
data follows Ge-ph ^ behavior (P cx T^) and fits well 
to the theory, but above 300 mK there is some deviation 
from the expected behavior. The origin of this devia- 
tion is not clear, but it is worthwhile noting that Eqs. (|3| 
are asymptotic zero Kelvin expressions for bulk phonon 
spectrum. 

The strained case {Fs) has no free parameters and 
we have plotted Ge-ph = dFs/dTe with dimensional- 
ity d = 2 and d = 3. These theoretical predictions fall 
significantly below our experimental data and in order to 
fit everything into same axis we have scaled dFs / dTe by 
a factor of 300. Here d = 2 means a quasi 2D-electron 
system where density response and screening are essen- 
tially 2D due to large thermal phonon wave length (that 
exceeds the electron layer thickness t), which we expect 
to be a better description of the system at hand. In this 
case the quasi 2D density of states U2 is approximated by 
^2 = tus^ where us is the 3D density of states. From Fig. 
[2|b) we observe that the theory is at least two orders 
of magnitude below the experimental data irrespective 
of the choice of d. In other words, the theory suggests 
that by introducing the strain we should have a factor 
of several 1000 decrease in Ge-ph ^ but we observe only 
a factor of 20 — 50 decrease in the temperature range 
T = 200 — 480 mK. On the other hand, the experimental 
Ge-ph of the strained sample follows the slope of the 2D 
theoretical curve {Ge-ph ^T^) 1 which suggests that by 
introducing a factor of ~ 30 larger deformation poten- 
tial constant (s) we would be able to explain the behavior 
of the strained sample. This would mean that the de- 
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formation potential constants should depend strongly on 
the strain, because the control sample fits to the theory 
so well (below 300 mK). We are not aware of strain af- 
fecting the deformation potential constants, but our data 
suggest the e-ph coupling mechanisms in heavily doped 
semiconductors requires further investigation. 

In conclusion, we have shown that by inducing strain 
to a degenerately doped Si layer, we have reduced its 
electron-phonon coupling by a factor of 20 — 50 in the 
temperature range 200 — 480 mK. Our results are in qual- 
itative agreement with theorjEI, but quantitatively the 
effect is smaller than expected. The reduced electron- 
phonon coupling in the strained Si might be beneficial 
for bolometric sensors and cooling applications. 
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